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Surface energies for different GaAs surface orientations have 
been calculated as a function of the chemical potential. We 
use an energy density formalism within the first-principles 
pseudopotential density-functional approach. The equilib- 
rium crystal shape (ECS) has been derived from the surface 
energies for the (110), (100), (111), and (111) orientations. 
Under As-rich conditions all four considered surface orienta- 
tions exist in thermodynamic equilibrium, in agreement with 
experimental observations. Moreover, our calculations allow 
us to decide on previous contradictory theoretical values for 
the surface energies of the (111) and (III) facets. 



I. INTRODUCTION 

The equilibrium crystal shape (ECS) is that shape 
which, in the limit of infinitely large volume, yields the 
minimum free energy of a crystal. For a given arbitrary 
surface orientation and unit cell the atomic reconstruc- 
tion that yields the lowest surface free energy can be 
determined. However, it is well known that in general 
this will not result in a thermodynamically stable situ- 
ation, because the surface can further lower its energy 
by faceting on a macroscopic scale. The ECS provides a 
set of surface orientations that exist in thermodynamic 
equilibrium. Except for some situations with degenerate 
surface energies surfaces of any other orientations will 
facet. 

The faceting of GaAs surfaces has been studied experi- 
mentally. Whereas Weiss et alB studied the different sur- 
face orientations exposed on a round shaped crystal jyith 
low-energy electron-diffraction (LEED) , Notzel et aln in- 
vestigated various planar high-index surfaces with re- 
flection high-energy electron-diffraction (RHEED). Both 
groups observed for different high-index surface orien- 
tations faceting into low-index surfaces. Moreover, sur- 
face energies play a major role in the formation of is- 
lands during heteroepitaxy. For example, InAs grows on 
GaAs in the Stranski-Krastanov mode.H The surface en- 
ergy of InAs being lower than that of GaAs, first a uni- 
form wetting-layer forms. During further deposition of 
InAs three-dimensional islands are formed due to strain 
relaxation. Receotly, these quantum dots have attracted 
great interest .Bu'ErQ Besides other quantities like the elas- 
tic relaxation energy of the islands, the absolute InAs 
surface energies of the involved facets which we assume 
to be similar to those of GaAs enter into the theory of 
the shape and size of the islands. 



Both experimental as well as calculated absolute val- 
ues of the surface energy as a function of orientation 
are quite scarce. The surface energy has been measured 
for the GaAs (110) surface in a fracture experiment.!!! 
Relative surface energies and the ECS of Si have been 
determinedou, but to our knowledge no such measure- 
ments have been carried through for GaAs. Moreover, 
it is often difficult to establish whether an observed sur- 
face really represents thermodynamic equilibrium. At 
low temperatures faceting and therefore thermodynamic 
equilibration may be hindered by insufficient material 
transport. At high temperatures, kinetics may govern 
the surface morphologies due to evaporation. 

The purpose of this work is to present the absolute val- 
ues for the surface energy of the GaAs (110), (100), (111), 
and (111) surfaces calculated from first principles, and 
the ECS constructed from these data. Empirical poten- 
tials do not produce reliable surface properties. Ab initio 
calculations have been carried out by various groups for 
different surface orientations of GaAs. Qian et al£3 used 
an ab initio pseudopotential method to calculate the ab- 
solute surface energy of the GaAs (110) surface. They 
found very good agreement with the. experimentaLcleav- 
age epergy. Northrup and FroyenEH, Qian et alH3, and 
OhnoEj determined the (100) reconstruction with low- 
est energy. The absolute surface energies for these 
constructions were not given, however. Kaxiras et al.t3 
calculated energies for GaAs (111) reconstructions rela- 
tive to the surface energy of the-ideal (111) surface. Foe 
the (111) surface Kaxiras et al.lB and Northrup et alJtB 
calculated relative surface energies for different (2 x 2) 
reconstructions. Based on their results they predicted 
the (111) equilibrium reconstruction. 

However, for geometrical reasons it is impossible to 
derive absolute surface energies for the (111) and (111) 



orientations of Ga As. from such total-energy calculations. 
Chetty and MartinE3'E£l solved this problem by introduc- 
ing an energy density, which enables the computation of 
the energies of the top and the bottom surfaces of the 
slab separately. Having calculated the absolute surface 
energies for the ideal reference surfaces they transformed 
the relative surface energies of Kaxiras et alE-$3 and 
Northrup et alx3 to absolute surface energies. A com- 
parison of these absolute values, however, shows that the 
two results differ significantly. This difference is not yet 
understood, and we will come back to it in Section IV 
below. 

We have calculated absolute surface energies for the 
different orientations directly (i.e., without introducing a 
reference surface) and consistently with one and the same 
set of parameters and pseudo-potentials. Before we will 
detail our results and the ECS of GaAs in Section IV, 
we will first give an overview of GaAs surface properties 
in Section II and describe the computational details in 
Section III. 



II. CHEMICAL POTENTIAL AND SURFACE 
RECONSTRUCTION 

The stable surface reconstruction is the one with the 
lowest surface free energy. In our case the substrate 
consists of two elements and thus the difference of the 
number of atoms of the two species enters as another 
degree of freedom in addition to the atomic geometry. 
Non-stoichiometric surfaces are considered by allowing 
the surface to exchange atoms with a reservoir, which is 
characterized by a chemical potential. The equilibrium 
is determined by the minimum of the free energy 

Tsurfacc^ — -^surface ^ ^ f^i^i- (1) 

i 

The surface free energy 7 su rfacc^4 of the surface area A 
has been calculated for zero temperature and pressure 
and neglecting zero point vibrations. The chemical po- 
tential fii is the free energy per particle in the reservoir 
for the species i, and Ni denotes the number of particles 
of the species i. The temperature dependence is ignored 
because the contributions tend to cancel for free energy 
differences. 

In experiment the value of the chemical potential can 
be varied over a certain interval. This interval is lim- 
ited by the bulk chemical potentials of the condensed 
phases of Ga and AsE3li3, corresponding to the two fol- 
lowing situations: On the one hand the surface can be in 
equilibrium with excess Ga-metal, which has the chemi- 
cal potential /^Ga(buik) j & n d the GaAs bulk with chemical 
potential /XGaAs- On the other hand the surface can be 
in equilibrium with bulk As and, again, the GaAs bulk. 
Both reservoirs can act as sinks and sources of surface 
atoms. The upper limit of each chemical potential is 



determined by the condensed phase of the respective el- 
ement, 

Hi < Mi(bulk): (2) 

because otherwise the elemental phase would form on the 
GaAs surface. Furthermore, in thermodynamic equilib- 
rium the sum of chemical potentials of Ga and As must 
be equal to the bulk energy per GaAs pair, 

^Ga + A* As = MGaAs 

= MGa(bulk) + A 1 As(bulk) ~ AHf. (3) 

For the heat of formation we have calculated a value of 

0. 64 eV using a plane-wave cutoff of 10 Ry_.which is in 
good agreement with the experimental valued of 0.74 eV. 
For the bulk calculations we computed the bulk energy of 
Ga in an orthorhombic structures and the bulk energy 
of As in a trigonal structurecd. 

In this work we give the surface energies in dependence 
of the As chemical potential. Therefore, we write equa- 
tions (||) and (||) in the following form 

MAs(bulk) - AHf < flAs < A 1 As(bulk) • (4) 

The surface energy is calculated from the total energy 
Etot, 

Tsurface^ = E tot - ^GaAs^Ga ~ MAs(A As - N Ga ). (5) 

The stoichiometry of the surface, AN — Na s — Nq &1 de- 
termines the slope of the surface energy versus the chem- 
ical potential. A consistent counting method for AN has 
to be applied to all orientations. We apply the method 
of Chetty and MartinEj which utilizes the bulk symme- 
tries of the crystal. For example, following their counting 
method the ideal (110) cleavage surface is stoichiometric, 

1. e. the difference AN is equal to zero. Thus the surface 
energy of the (110) cleavage surface is independent of the 
chemical potential. 

When the chemical potential is varied, different recon- 
structions with different surface stoichiometrics become 
thermodynamically stable. All experimentally observed 
reconstructions, however, fulfill certain conditions. First 
of all, GaAs surfaces favor to be semiconducting, as this 
leads to a low surface energy. Surface bands in the bulk 
gap and especially surface bands crossing the Fermi-level 
will lead to,a-higher surface energy. The electron count- 
ing modelE3a gives a simple criterion whether a surface 
can be semiconducting or not. In the bulk the sp 3 hy- 
bridized orbitals of GaAs form bonding and antibonding 
states. At the surface there are partially filled dangling 
bonds. Their energies are shown schematically in Fig. |lj, 
they are estimated from the atomic s and p eigenener- 
gies of either species. Compared to the dispersion of the 
conduction and the valence bands, the dangling bond en- 
ergy of the cation (Ga) falls into the conduction band and 
therefore it should be empty. The dangling bond energy 
of the anion (As) lies in the valence band and therefore 
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FIG. 1. The energy levels of the s and p orbitals e s , p , of 
the sp 3 dangling-bonds eh, and of the conduction and valence 
band. The data are from Harrison.E3 



it should be filled. Thus there has to occur an electron 
transfer from the Ga to the As dangling bonds. For a 
low-energy semiconducting surface the dangling bonds 
in the conduction band have to be empty, exactly filling 
all the dangling bonds in the valence band. Otherwise 
the surface becomes metallic and has a higher surface 
energy. Ga and As surface atoms are added to, or re- 
moved from, the ideal bulk-truncated polar surfaces to 
obtain a low-energy semiconducting surface. 

Secondly, the electron transfer from the Ga dangling 
bonds to the As dangling bonds has consequences for 
the geometry of the surface reconstructions. The surface 
Ga atom which has lost an electron favors a sp 2 like hy- 
bridization. Therefore the Ga atom relaxes inwards and 
forms a more planar configuration. The dangling bond 
of arsenic is completely filled and the As atom prefers 
to form bonds with its three p orbitals. Therefore the 
bond angle of the surface As atom is close to 90°, and 
the As atom relaxes outwards. These configurations re- 
semble the bond geometry of small molecules like GaH 3 
and ASH3 and axe a general result for surfaces of III-V 
semiconductors O 



III. COMPUTATIONAL DETAILS 

To determine the surface energies we carried out tptaA-, 
energy calculations using density-functional theoryESES 
We applied the local-density approximation to the 
exchange-correlation functional, choosing the parame- 
terization by Perdew and ZungercH of Ceperley and 
Alder'sca data for the correlation energy of the homo- 
geneous electron gas. The surfaces were described by 
periodically repeated slabs. All computations were done, 
with an extended version of the computer code fhi93cpE21 



The program employs ab initio pseudopotentials and a 
plane- wave basis-set. It was generalized to addition- 
ally cornnute the energy density according to Chetty and 
MartinO 

The slab geometry leads to serious problems when sur- 
face energies of zinc-blende structures are to be calcu- 
lated for arbitrary orientations. To derive the surface 
energy from a total energy calculation both surfaces of 
the slab have to be equivalent. Though such slabs can be 
constructed for the (110) and the (100) orientation, this 
is impossible for the (111) orientation: The (111) and 
the (111) surfaces of GaAs arc inequivalcnt. This follows 
from the simple geometric property of the zinc-blcndc 
structure that the Ga-As double layers are Ga and As 
terminated on the top and bottom side of the slab, re- 
spectively. Chetty and MartirO solved this problem by 
introducing an energy density. The energy density itself, 
however, does not bear any physical significance, only the 
integrals of the energy density over suitable parts of the 
supercell (e.g., volumes bounded by bulk mirror plaxics) 
lead to well-defined, physically meaningful energies. t3 We 
have checked the accuracy of this approach for our GaAs 
slabs: Variation of the surface reconstruction on the bot- 
tom side of the (100) and (111) slabs results in a negli- 
gible change of the surface energy of the surface on the 
top (< 0.7 meV/A 2 ). 

Ab initio norm-conserving pseudopotentials were gen- 
erated with Hamann's schemc.Ej The cutoff radii for 
pseudoization have been chosen equal to 0.58 A, 0.77 A, 
and 1.16 A for the s, p, and d wave-functions of Ga, 
and equal to 0.61 A, 0.60 A, and 1.07 A for s, p, and 
d wave-functions of As. The semi-local pseudopotentials 
were further transformed into fully separable Klcinman- 
Bylander pseudopotentialsEl, with the d potential cho- 
sen as the local potential. The logarithmic derivatives of 
the different potentials were examined and various trans- 
ferability testsE 2 ], e.g. "hardness" tests, were performed. 
All together the potentials showed good transferability. 
The structures of the bulk phases of Ga and As are well 
described by these potentials, the theoretical lattice con- 
stants being only slightly smaller than the experimental 
ones with a relative deviation below 3.5%. ._. 

The wave functions were expanded into plane wavesEj 
with a kinetic energy up to 10 Ry. This leads to a conver- 
gence error in the surface energies of less than 3 me v/A 2 . 
The-electron density was calculated from special k-point 
setsc-a, their density in reciprocal space being equivalent 
to 64 k-points in the whole (100) (1 x 1) surface Brillouin- 
zone. 

For the (100), (111), and (111) surfaces "pseudo- 
hydrogen" was used to saturate the bottom surfaces 
of the slabs £3 Pseudo- hydrogen denotes a Coulomb- 
potential with a non-integer core-charge Z, together with 
Z electrons. The Ga and As atoms of these surfaces 
were fixed at their ideal bulk positions. The Ga termi- 
nated surface was saturated with pseudo-hydrogen with 
an atomic number of Z — 1.25. On each dangling bond 
of a Ga surface atom one pseudo-hydrogen was placed. 
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Similarly, the As terminated surface was saturated with 
pseudo-hydrogen with an atomic number of 0.75. The 
saturated surfaces are semiconducting without any sur- 
face states in the bulk band-gap. There are two main 
advantages using this pseudo- hydrogen. First of all, the 
interaction of both surfaces with each other is in this 
way minimal. Secondly, the surface atoms which are sat- 
urated with the pseudo-hydrogen can be kept fixed at 
ideal bulk positions. Thus thinner slabs can be used and 
charge sloshing is suppressed. 

For polar surfaces, such as the ideal (111) surface, a dif- 
ficulty arises due to charge transfer from one side of the 
slab to the opposite side. This charge transfer is hindered 
by a semiconducting surface, e.g. the pseudo-hydrogen 
saturated surface at the bottom of the slab. We esti- 
mate the uncertainty due to charge transfer to be smaller 
than 1.4 meV/A 2 for a polar surface, comparing the sur- 
face energies of the pseudo-hydrogen saturated surface 
derived from two calculations. One is carried through 
with a semiconducting surface on the top of the slab, the 
other one with a metallic surface. 

We have carried out computations for a large vari- 
ety of reconstructions of the GaAs (110), (100), (111), 
and (111) surfaces, which have previously been suggested 
in literature. Starting from some initial geometry, the 
atom positions in the topmost layers of the slab were re- 
laxed until the forces on the atoms were smaller than 50 
meV/A. The other layers were kept fixed at their ideal 
bulk positions with a bulk lattice-constant of 5.56 Awhich 
had been determined theoretically at the same cutoff en- 
ergy as the slab calculations and using 384 k-points in 
the whole Brillouin-zone. This value is 1.4% smaller than 
the experimental lattice constants neglecting zero point 
vibrations. 



IV. RESULTS AND DISCUSSION 

A. (110) Surface 

The (110) surface is one of the most extensively studied 
GaAs surfaces (Ref. [n| |37|, [3g, pijand references therein). 
The (110) plane is the cleavage plane of III-V semicon- 
ductors. Containing the same number of cations (Ga) 
and anions (As) it is intrinsically neutral. The cleav- 
age surface does not reconstruct, only a relaxation of 
surface atomic positions within the (lxl) surface unit 
cell is observed. The charge from the Ga dangling-bond 
is transfered into the As dangling-bond, which becomes 
completely filled. The orbitals of both surface atoms re- 
hybridize, and the zigzag chains of Ga and As surface 
atoms tilt, with the As atom being raised and the Ga 
atom being lowered. Thereby the Ga surface atom ac- 
quires a nearly planar bonding configuration, while the 
As surface atom relaxes towards a pyramidal configura- 
tion with orthogonal bonds. 




(a) cleavage 




(b) Ga terminated (c) As terminated 



FIG. 2. Atomic structures of the GaAs (110) surface in 
top and side view. Open and filled circles denote As and Ga 
atoms, respectively. 

We have calculated the surface energy of the relaxed 
cleavage surface shown in Fig. ||(a). It is stoichiometric 
(AN = 0) and semiconducting. In addition, we consid- 
ered two other surface structures: The Ga terminated 
(110) surface is shown in Fig. ||(b). Formally it can be 
constructed from the cleavage surface by substituting all 
top-layer As atoms by Ga atoms. This surface has a sto- 
ichiometry of AA" = —2 per (1 x 1) cell, and it fulfills 
the electron counting criterion. Nevertheless, it is not 
semiconducting, because the bands of the Ga-Ga surface 
bonds cross the Fermi-level. The Ga surface atoms do 
not relax in the same way as the respective Ga and As 
atoms in the cleavage surface, instead they almost stay 
in the same plane. Finally, we have calculated the sur- 
face energy of the As terminated (110) surface (see Fig. 
H(c)). Here the Ga surface atoms have been replaced 
by As atoms, which yields a surface with a stoichiome- 
try of A A" = 2 per (lxl) surface unit cell. Also this 
surface fulfills the electron counting criterion, and it is 
semiconducting. Both As dangling-bonds are completely 
filled and lie beneath the Fermi-level. Similar to the Ga 
terminated surface, the As surface atoms do not relax 
significantly, but stay in the same plane. 

For all three (110) surface reconstructions we used the 
same super cell, with slabs composed of nine atomic lay- 
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M-As M-As(bulk) ( e ^) 

FIG. 3. Surface energy of the different GaAs (110) sur- 
face structures in me V/A 2 plotted versus the difference of 
the chemical potential of As and As bulk. 

ers and a vacuum region with a thickness equivalent to 
seven atomic layers. The whole surface Brillouin zone 
was sampled with 48 special k-points.E.3 

The calculated surface energies are shown in Fig. || for 
the three surface structures we have considered. For a 
large range of the chemical potential the cleavage surface 
is energetically most favorable. Our result for the sur- 
face energy of 52 meV/A 2 is in good agreement with the 
value of 57eV/A 2 meV which was calculated by Qian et 
alS-i using essentially the same ab initio method. Both 
results compare very well with the experimental surface 
energy of 54 ± 9 meYj/A 2 from fracture experiments by 
Messmer and Bilello.0 In As-rich environments we find 
the As terminated surface to exist in thermodynamic^], 
equilibrium, in agreement with Northrup's calculation.Ej 
We obtain a value of 45 meV/A 2 for the. surface energy in 
an As-rich environment. Kiibler et alEU provided experi- 
mental evidence for the existence of this structure. Using 
LEED they observed that the surface relaxation was re- 
moved as the As coverage was increased. In contrast to 
the As terminated surface, we find the Ga terminated 
surface to be unstable even under the most extreme Ga 
rich conditions. 



B. (100) Surface 

Among the different orientations the (100) surface 
is the one used most widely for the growth of opto- 
electronic devices. The (100) surface is polar, i.e. the 
planes parallel to the surface consist of either only Ga 
or only As atoms. As a consequence, the stable sur- 
face structures displays various reconstructions which 
distinctly differ from those found on the (100) faces p£ 
the covalent group IV semiconductors. Daweritz et aZ.E£l 



have derived a steady state "phase" diagram for the sur- 
face reconstruction as a function of growth conditions. 
In their diagram they point out 14 different reconstruc- 
tions. To our knowledge, the equilibrium phase diagram 
of the (100) surface has not yet been determined. How- 
ever, there are certain reconstructions which are gener- 
ally observed during anxLafter growth. While heating the 
surface Biegelsen et a/Jl2l observed a sequence of phases 
from the As-rich c(4 x 4), (2x4) to the Ga-rich (4 x 2) 
reconstructions. For each of these surface unit-cells there 
exists a isjrge variety of possible atomic configurations. 

Chadicil performed tight binding based total energy 
minimizations to examine the structure of the (2x1) 
and (2 x 4) reconstructed surface. For the (2 x 4) he sug- 
gested two possible atomic configurations with three and 
two As.-dimers (f3 and (32, notation according to Northrup 
et alEn) per surface unit-cell. Moreover, he determined 
the energy difference between the (2 x 4) and the related 
c(2 x 8) reconstruction to be less than 1 meV/A 2 . As 
the (2 x 4) and the c(2 x 8) are very similar and have 
only small difference in surface energy, we have not cal- 
culated, the centered reconstructions c(2 x 8) and c(8 x 2). 
Ohno£3 and NorthrupO carried through ab initio calcula- 
tions of the surface energies. Ohno could exclude various 
configurations of the (2x1) and (3x1) surface unit-cell 
from being equilibrium structures. Moreover, he con- 
cluded that for the (2 x 4) reconstruction the phase j3 
with three surface dimers is stable, which appeared to be 
in-agreement with the STM observations of Biegelsen et 
a/o However, calculations by Northrup et a/.Ej showed 
that the most stable (2 x 4) reconstruction contains two 
As-dimers in the top layer, which has beep, confirmed 
by recent high resolution STM observations E3 Northrup 
et al. also investigated the energetics of the (4 x 2) and 
c(4 x 4) reconstructions. For the (4 x 2) reconstruction 
they found a two-dimer phase to be energetically favor- 
able in agreement with STM investigations O Howevec- 
a recent analysis of LEED intensities by Cerda et als3 
suggests that the top layer consist of three Ga dimers 
per (4 x 2) unit cell. For the c(4 x 4) reconstruction 
Northrup et al. considered a three-dimer phaseo which 
they found to be stable in certain conditions with re- 
spect to the (2 x 4) and (4 x 2) reconstructions. On the 
other hand a ±wo-dimer phase was suggested by Sauvage- 
Simkin et a/.cfj on the basis m£ X-ray scattering experi- 
ments, and by Larsen et alrA who studied the surface 
with a number of different experimental techniques. 

In our calculations we have considered all atomic con- 
figurations with a (2 x 4) and a (4 x 2) surface unit-cell, 
that were previously investigated by Northrup et a/O'C 2 ] 
For the c(4 x 4) reconstruction we took into account the 
three-dimer phaseo and a structure which has two in- 
stead of three As-dimers in the top layer.E3c3 The to- 
tal energy calculations were performed using supercells 
containing seven layers of GaAs. The thickness of the 
vacuum corresponded to five layers GaAs. 

In Fig. U the geometries of those surface structures are 
shown that have minimum surface energy within some 
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FIG. 4. Atomic structures of the GaAs (100) reconstruc- 
tions. 



range of the chemical potential and therefore exist in 
thermodynamic equilibrium. All four structures fulfill 
the electron counting criterion and are semiconducting, 
i.e., the anion dangling bonds are filled and the cation 
dangling bonds are empty. Furthermore the surfaces dis- 
play Ga-Ga bonds and As- As bonds, both having filled 
bonding and empty antibonding states. The a (2 x 4) re- 
construction (Fig. ||(a)) is stoichiometric (AN = 0). In 
the top layer four As atoms are missing per (2 x 4) cell. 
The surface As atoms form two dimers. The Ga-layer un- 
derneath is complete, but differs from the bulk geometry 
by two Ga-Ga bonds which are formed between the Ga 
atoms in the region of the missing As dimers. Removing 
the Ga atoms in the missing dimer region one obtains 
the /32(2 x 4) structure in Fig. f|(b) with a stoichiome- 
try of AA~ = j per (lxl) unit cell. The completely 
As-terminated c(4 x 4) surface shown in Fig. ||(c) has a 
stoichiometry of AA?" = | per (1 x 1) unit cell. It consists 
of three As-dimers which are bonded to a complete As- 
layer beneath. The /32(4x 2) structure shown in Fig. [|(d) 
represents the Ga-terminated counterpart of the /32(2x4) 
reconstruction, with Ga atoms exchanged for As atoms 
and vice versa. Thus the top layer consists of two Ga 
dimers per (4 x 2) cell, and the second layer lacks two As 
atoms. This results in a stoichiometry of AA" = —j per 
(1 x 1) cell. 

Our calculated surface energies of these four phases 
are shown in Fig. (^ as a function of the chemical poten- 
tial. We predict the same sequence of equilibrium surface 
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M-As - M-As(bulk) ( e ^) 

FIG. 5. Surface energy of the different GaAs (100) recon- 
structions in meV/A 2 plotted versus the difference of the 
chemical potential of As and As bulk. 



structures as Northrup and Froyent-UO as a function of 
increasing As coverage: /32(4x2), a(2x4), /32(2x4), and 
c(4x4). The c(4 x 4) structure with only two surface 
As-dimers per unit cell, which we considered in addition 
to the structures investigated by Northrup and Froyen, 
turned out to be unstable. Though this structure is more 
Ga-rich than the c(4x4) three As-dimer structure shown 
in Fig. ^(c), even in the Ga-rich environment the two- 
dimer phase has a surface energy which is 5 meV/A 2 
higher than for the three-dimer phase. Due to the lack 
of absolute values in previous calculations, quantitatively 
we can only compare energy differences between surfaces 
with the same stoichiometry. Further comparison is made 
difficult by the different range of the chemical potential in 
our versus Northrup and Froyen's calculatiorHo: Their 
value for the heat of formation is AHf = 0.92 eV, as op- 
posed to our smaller value of AHf = 0.64 eV. Compar- 
ing the three dimer phase (3 with the two dimer phase /32, 
which both have the same stoichiometry, we find that the 
two dimer phase has a surface energy lower by 2 meV/A 2 . 
This agrees with the result of Northrup and Froyen who 
report an energy difference of 3 meV/A 2 , and it further 
confirms the conclusion that the three dimer phase (3 does 
not exist in equilibrium. On the whole, the agreement 
with the relative surface energies calculated by Northrup 
et al. is good. They can be converted to absolute surface 
energies by shifting them by w 65 meV/A 2 , which results 
in a diagram similar to Fig. ||. 

All investigated (100) surfaces display similar atomic 
relaxations which are characterized by the creation of 
dimers and the rehybridization of threefold coordinated 
surface atoms. The creation of surface dimers decreases 
the number of partially occupied dangling bonds, and by 
rehybridization the surface gains band structure energy. 
The calculated bond lengths in bulk Ga and As, 2.32 A 
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and 2.50 A, respectively, can serve as a first estimate for 
the respective dimer bond lengths on the GaAs surface. 
Our calculations yield As-dimer lengths between 2.45 and 
2.50 A for the a and 02 surface reconstruction. This is 
within the range of experiment all y. .deduced values which 
scatter between 2.2 and 2.9 AoO'Ej and it is similar to 
the dimer lengths of 2.53-and 2.55 A which were deter- 
mined by Northrup et alE3 On the c(4 x 4) reconstructed 
surface the calculated As-dimer lengths are 2.57 A for 
the central dimer and 2.53 A for the two outer dimers 
of the three-dimer strings in the surface .writ-cell. Us- 
ing X-ray scattering Sauvage-Simkin et alc2 determined 
these bond-lengths as 2^3 ± 0.06 A and 2.59 ± 0.06 A. 
Very recently Xu et alx2 suggested that the dimers on 
the c(4x4) structure should be tilted by 4.3°. However, 
as for the (2 x 4) reconstructions we find the dimers to be 
parallel to the s urfa ce, in agreement with several previ- 
ous experimentsc3c3. Even when starting with an initial 
configuration with surface dimers tilted by 8° we find 
the dimers to relax back to the symmetric positions with 
a residual tilt angle less than 0.1°. The Ga-Ga dimer 
bond length is calculated to be 2.4 A on the /32(4 x 2) 
reconstruction and 2.5 A on the a(2 x 4) structure which 
agrees with previous ab initio calculations .EiJ From a re- 
cent LEED investigation of the Ga rich (100) surface 
Cerda et alnB deduced that the stable (4 x 2) recon- 
structed surface displays three Ga-dimers per unit cell 
with unusual dimer lengths of 2.13 A and 3.45 A. In our 
calculation, however, this three dimer phase is energet- 
ically slightly less favorable than the two dimer phase 
/32(4 x 2) by 0.8 meV/A 2 . Therefore, it should not be 
stable at least at low temperatures. Furthermore, we 
found the Ga dimer length to be 2.4 A and no local min- 
imum for Cerda's unusually large dimer length. 

The rehybridization of the sp 3 orbitals located at the 
threefold coordinated Ga-atoms drives the relaxation to- 
wards a preferentially flat Ga-bond configuration. On 
the Ga terminated /?2(4 x 2) structure this leads to a de- 
creased spacing between the Ga top layer and the neigh- 
boring As layer which amounts to roughly half of the bulk 
interlayer spacing. Also on the a(2 x 4) and /32(2 x 4) sur- 
faces the threefold coordinated Ga atoms which bond to 
As relax towards the plane of their neighboring As atoms. 
Together with a slight upward shift of the top layer As 
atoms this leads to a steepeningE3 of the As dimer block. 
The change of the angle between the bonds of the three- 
fold coordinated As atoms is less pronounced. However, 
the trend is obvious: except for the c(4 x 4) structure, 
we find the As bond-angles to be always smaller than 
109.5°, which is the angle of the ideal tetrahedral co- 
ordination. The As bonds on the c(4 x 4) surface be- 
have differently from those on the other three surfaces 
because the top-layer As atoms are bonded to a second 
layer which consist of As instead of Ga. A decrease of the 
angle between the bonds of all threefold coordinated As 
atoms would require a change in the As- As bond lengths, 
which probably costs more energy than would be gained 
from rehybridization. 




(a) ideal (b) As adatom 




(c) As trimer (d) Ga vacancy 

FIG. 6. Atomic structures of the GaAs (111) reconstruc- 
tions. 

C. (Ill) Surface 

The polar (111) orientation of GaAs has been .slucLied 
within density-functional theory by Kaxiras et a/HoEf!, 
who computed surface energies relative to the ideal unre- 
constructed surface for various atomic geometries. They 
found that under As-rich conditions an As trimer ge- 
ometry yields the lowest surface energy, whereas a Ga 
vacancy reconstruction is prefaced under Ga-rich condi-, 
tions. Haberern and PashleyEj and Thornton et alxQ 
confirmed this experimentally. Haberern and Pashlcy 
interpreted their STM images to show an array of Ga 
vacancies with a (2x2) periodicity. Thornton et al. ob- 
served both the As triangle model and the Ga vacancy 
model in STM. Here we concentrate on the following re- 
constructions of the Ga terminated (111) surface: the As 
adatom, the As trimer, the Ga vacancy model, and, for 
comparison but not as a reference system as in previous 
work, the truncated-bulk geometry. 

The ideal (111) surface (see Fig. ^(a)) has a stoichiom- 
etry AN = — 4/(1 x 1). It does not fulfill the electron 
counting criterion. Each Ga dangling-bond is filled with 
3/4 of an electron and therefore the ideal surface has to 
be metallic. To create a neutral semiconducting surface, 
following the electron counting criterion one can either 
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add an As surface atom to, or remove a Ga surface atom 
from, every (2x2) surface unit cell. Therefore we consider 
three different (2 x 2) reconstructions. First of all, the As 
adatom model is shown in Fig. |^(b) . This reconstruction 
is stoichiometric. The As adatom binds to the Ga surface 
atoms. It exhibits a nearly orthogonal bond configura- 
tion, while the Ga atom with the empty dangling bond 
relaxes towards the plane of the As atoms. Secondly, we 
consider the As trimer model shown in Fig. |6|(c). This 
model has a stoichiometry AN = i/(lx 1), it also fulfills 
the electron counting criterion and it is semiconducting. 
The three extra As atoms form a trimer with each As 
atom binding to one Ga atom. The dangling-bonds of 
the As atoms are completely filled and the dangling-bond 
of the the Ga atom which is not bonded to As trimer 
atoms is completely empty. This Ga atom relaxes into 
the plane of the As atoms of the layer below. Finally, 
we calculated the Ga vacancy model (see Fig. ||(d)). The 
removal of one Ga surface atom causes the surface to be 
stoichiometric. The Ga surface atoms have completely 
empty dangling-bonds and relax into the plane of the 
As atoms. The three As atoms surrounding the vacancy 
have completely filled dangling-bonds. 

We used the same super cell for the calculations of the 
(111) and the (111) surfaces. Only the bulk-truncated 
surface was calculated within a (1 x 1) surface unit cell, 
else always a (2 x 2) unit cell was used. The slab con- 
sisted of five (111) double layers. The vacuum region had 
a thickness equivalent to four (111) double layers. The 
whole Brillouin zone of the (2 x 2) surface unit cell was 
sampled with 16 special k-points, corresponding to 64 k- 
points in the Brillouin zone of the (1 x 1) cell. Absolute 
surface energies of the (111) reconstructions were deter- 
mined using the energy density formalism. The results 
are shown in Fig. [t]. The Ga vacancy model is the most 



favorable reconstruction for a large range of the chemi- 
cal potential from a Ga-rich to an As-rich environment. 
Only in very As-rich environments the As trimer model 
has a lower energy. The Ga vacancy model has a surface 
energy of 54 meV/A 2 , whereas, the As trimer model has 
a surface energy of 51 meV/A 2 in As-rich environment 
at //As — / i As(buik) ■ The Ga vacancy reconstruction was 
observed experimentally by Haberern and PashleyE3 and 
Tong et alxli Thornton et al. additionally observed the 
As trimer reconstruction. 

Two other groups have performed similar ab initio cal- 
culations. Using their energy density formalism, Chetty 
and Martinlia derived a value of 131 me V/A 2 for the sur- 
face energy of the ideal (111) surface in a Ga-rich en- 
vironment, which is much larger than our value of 93 
nic V/A 2 . Secondly, we can compare outresults to the rel- 
ative surface energies of Kaxiras et alE3^3 They arrived 
at the same qualitative conclusions. However, quantita- 
tively their relative surface energies are not easily com- 
parable to ours because they used AS4 gas to define the 
As-rich environment. Therefore they obtained a larger 
interval for the chemical potential. We derive for the sur- 
face energy difference of the As adatom and Ga vacancy 
structure a value of 13 meV/A 2 , whereas Kaxiras et al. 
calculate a much larger difference of 47 me V/A 2 . Using 
their own result for the ideal surface, Chetty and Martin 
transformed the relative surface energies of Kaxiras et al. 
to absolute surface energies. In comparison to our results, 
all these surface energies contain the same shift towards 
higher energy as the ideal surface mentioned above. We 
will discuss this difference below and explain, why we 
believe our results to be accurate. 

Tong et alEli performed a LEED analysis for the ge- 
ometry of the Ga vacancy reconstruction. Their geom- 
etry data_.compare very well with the theoretical data 
of ChadiEj, Kaxiras et aZo and ours. For the Ga va- 
cancy reconstruction we find an average bond angle of 
the sp 2 -bonded Ga surface atom of 119.8° in agreement 
with Tong et al. The bond angles of the p 3 -bonded As 
atom of 87.0° and 100.6° average to 91.5 ° which again 
compare very well with the value of 92.9° by Tong et 
al. The bonds of the p 3 -bonded As atom are strained by 
-1.6% and 2.6% with respect to the GaAs bulk bonds. 
Tong et al. measured a value -1.3% and 1.9%, respec- 
tively. 

Furthermore, for the As trimer reconstruction we corn- 
pace our geometry data to theoretical data of Kaxiras et 
alE3 The threefold-coordinated As adatoms form bond 
angles to the neighboring As adatoms of 60° due to sym- 
metry reasons. The bond angle of the As adatom to 
the next Ga atom is 106.2°. Therefore we get an aver- 
age bond angle of 90.8° which is in good agreement with 
the 91.7° of Kaxiras et al. The surface Ga-As bonds are 
strained by 1.4 %, whereas Kaxiras et al. find the same 
bond length as in the bulk. The As- As bonds have a bond 
length of 2.44 A, 2.4 % shorter than that in As bulk. The 
Ga surface atom which is not bond to an As adatom re- 
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laxes into the plane of the As atoms with a bond angle 
of 118.4° and a bond length which is 2.6 % shorter than 
in GaAs bulk. These values are slightly larger than the 
114.7° and 1.0 % reported by Kaxiras et al. 

D. (Ill) Surface 

The polar GaAs (111) surface differs from the (111) 
surface, as the bulk-truncated (111) surface is terminated 
by As atoms, while the (111) surface is Ga terminated. 
At first sight the (111) surfaces might seem to be still 
analogous to the (111) surfaces, only that the Ga and As 
atoms have to be exchanged. However, this analogy is 
not useful, because As and Ga have different electronic 
properties, and therefore the (111) and (111) surfaces 
do not exhibit equivalent reconstructions. Stoichiometric 
(TIT) surfaces are gained by adding a Ga atom per (2 x 
2) surface unit cell to the bulk-truncated surface or by 
removing an As_surface atom. 

Kaxiras et aZji3 calculated the relative surface enerps 
for various (2 x 2) reconstructions. Biegelsen et al.tB 
studied the (111) surface both experimentally and theo- 
retically. Using STM they observed an As trimer (2x2) 
reconstruction for As-rich environments. A (\/l9 x \/l9) 
reconstruction which is dominated by two-layer hexago- 
nal rings was identified for Ga-rich environments. 

Due to the large unit cell the (\/l9 x VT9) reconstruc- 
tion is computationally quite expensive, and in this work 
we thus only consider (2 x 2) reconstructions. First of 
all, for comparison, we calculate the surface energy of the 
ideal (i.e., relaxed bulk-truncated) surface shown in Fig. 
||(a). This surface is not stoichiometric (AN — j/(lxl)). 
The dangling-bond of each As surface atom is filled with 
5/4 of an electron. Therefore the surface is metallic. 
Secondly, the Ga adatom model shown in Fig. ||(b) was 
considered. Through adding of an additional Ga surface 
atom the surface has become stoichiometric and semicon- 
ducting. The dangling-bond of the Ga adatom is com- 
pletely empty, whereas the dangling-bond of the As atom 
which is not bond to the Ga adatom is completely filled. 
Furthermore, we also consider an As trimer model (see 
Fig. H(c)). In contrast to the (111) surface the As trimer 
is bond to As surface atoms. This reconstruction has a 
stoichiometry of AN = 1 per (lxl) surface unit cell. 
Each As surface atom has a completely filled dangling- 
bond. Therefore, the surface is semiconducting. Further- 
more, we calculate the surface energy for the As vacancy 
model which is shown in Fig. ||(d). The removal of the 
As surface atom causes the surface to be stoichiometric. 
The three neighboring Ga atoms have completely empty 
dangling-bonds. The surface fulfills the electron counting 
criterion and is semiconducting. Finally, we calculate the 
Ga tipper model (see Fig. ||£&}) to compare with Kaxiras 
et alr3 and Northrup et alS3 This surface model has a 
stoichiometry of AN = — 1/2 per (lxl) surface unit cell 
and also fulfills the electron counting criterion. However, 




(a) ideal (b) Ga adatom 




(c) As trimer (d) As vacancy 




(e) Ga trimer 

FIG. 8. Atomic structures of the GaAs (Til) reconstruc- 
tions. 

it is metallic for the same reason as the Ga terminated 

(110) surface. 

The calculations for the (TIT) surface were carried out 
with the same parameters and supercell as those for the 

(111) surface outlined in the previous section. The re- 
sults are shown in Fig. ^|. For As-rich environments we 
find that the As trimer model is the most favorable re- 
construction, as observed experimentally by-STM and 
confirmed by previous ab initio calculations.^ This re- 
construction has a very low surface energy of 43 me v/A 2 . 
In a Ga-rich environment the Ga adatom reconstruction 
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has the lowest energy (69 meV/A 2 ) among all the struc- 
tures we calculated. The (\/T9 x \/l9) reconstruction 
found experimentally was not included in our approach. 
However, as suggested by Biegelsen et al. our present 
data can be used to restrict the range of possible values 
for the surface energy of the (VT9 x VT9) reconstruction 
consistent with observation: It has to be smaller than the 
surface-energy of the Ga adatom model on the one hand, 
and it has to be larger than the minimum energy of the 
As-trimer surface (plus a small correction of —3 meV/A 2 
to account for the non-stoichiometry of the \/l§x y/l§ re- 
construction) on the other hand. Therefore, we conclude 
that energy of the (VT9 x \/l9) reconstruction is in the 
range between 40 and 69 meV/A 2 . Considering also the 
energetical competition with facets of other orientations, 
even a slightly more stringent condition can be deduced: 
For the (111) (vl9 x \/l9) surface in a Ga-rich environ- 
ment to be stable against faceting into {110} surfaces, its 
surface energy has to be less than 63 meV/A 2 . 

In comparison to the relative surface energies cal- 
culated by Kaxiras et aZUj our energy difference be- 
tween the As vacancy and the Ga adatom structure of 
2 meV/A 2 is only slightly smaller than their value of 6 
meV/A 2 . However, they state that for the Ga-rich en- 
vironment the Ga trimer structure is 24 me v/A 2 more 
favorable than the Ga adatom structure. In contrast, we 
agree with Northrup et alS3 that the Ga trimer is ener- 
getically quite unfavorable. It has a 29 meV/A 2 higher 
surface energy than the Ga adatom. Also the other rela- 
tive surface energies compare quite well mith the already 
mentioned calculations of Northrup et alS3 although they 
derived a larger heat of formation (0.92 eV as opposed to 
our value of 0.64 eV). Relative to the Ga adatom our sur- 
face energies of the As trimer are about 10 me V/A 2 lar gcr 
than theirs. Also, they find a slightly larger energetic sep- 



aration between the As vacancy and Ga adatom struc- 
tures. Their value for this energy difference is 6 me v/A 2 , 
whereas our value is 2 meV/A 2 . However, these differ- 
ences are small and do not affect the physical conclusions. 
Chctty and Martin derived the absolute surface energies 
using their result for the ideal (111) surface and the rel- 
ative surface energies of Kaxiras et al. and Northrup et 
al. In contrast to the (111) their value of 69 meV/A 2 
for the ideal (111) surface in the Ga-rich environment is 
much smaller than ours of 97 mc V/A 2 . Therefore, this 
time in comparison to our data the results all contain the 
same shift to lower surface energies as the ideal surface. 
However, the sum of the (111) and (111) surface ener- 
gies from Chetty and Martin is close to ours. Therefore 
it is the splitting of the slab total energy into contribu- 
tions from the (111) and the (111) side that comes out 
differently. In our calculations both sides are energeti- 
cally similar which seems to be plausible in view of the 
fact that the flat (i.e., not faceted) surfaces have been 
observed experimentally. 

With respect to the calculated geometry we find that 
the As- As bond length in the trimer is 2.46 A, 1.6 % 
shorter than in bulk As. The As trimer atoms each bind 
to an As atom 2.30 A beneath, the As trimer plane in 
agreement with Northrup et al&3 The remaining As atom 
which is not bond to the trimer relaxes outwards and is 
1.74 A below the trimer plane. This compares reasonably 
well with the slightly larger value of 1.89 A by Northrup 
et al. For the two Ga surface models the separation of 
the adatom or trimer plane and the closest As (rest atom) 
plane amounts to 0.98 A for the Ga adatom model, and 
1.98 A for the Ga trimer model. Northrup et al. derived 
values of 0.98 A and 1.90 A. 



E. Equilibrium Crystal Shape (ECS) 

As opposed to liquids, crystals have non-trivial equi- 
librium shapes because the surface energy 7(m) depends 
on the orientation m of the surface relative to the crys- 
tallographic axes of the bulk. Once 7(111) is kapwn, the 
ECS is determined by the Wulff constructionists which 
is equivalent to solving 



r(h) 



7(m) 
™ Vih • h 



(6) 



Here r(h) denotes the radius of the crystal shape in the 
direction h. When the surface energy 7(1x1) is drawn as 
polar plot, the ECS is given by the interior envelope of 
the family of planes perpendicular to m passing through 
the ends of the vectors 7(111) rh. Under the assumption 
that only the (110), (100), (111), and (111) facets exist, 
we construct the ECS from the calculated surface en- 
ergies of these facets. Thus there may exist additional 
thermodynamically stable facets that are missing on our 
ECS. To be sure to construct the complete shape one 



10 



(010) 





w 


(110) 


(111) 


(001) 1 




1 




(101 
A 



(001) 



(100) 

FIG. 10. Three-dimensional representation of the ECS of 
GaAs in an As-rich environment, constructed from the surface 
energies of the (110), (100), (111), and (III) facets. The 
(l00), (010), and (001) axes are drawn for convenience. 



would have to calculate the surface energy for every ori- 
entation. However, from experiments it is known that 
the low Miller-indices surface orientations we consider 
are likely to be the energetically most favorable ones. 

As the GaAs surface energies depend on the chemical 
environment, the ECS becomes a function of the chemi- 
cal potential. In Fig. [l(] the ECS is shown for an As-rich 
environment and zero temperature. The different facets 
have been marked in the Figure and the ECS reflects 
the symmetry of bulk GaAs. To investigate the depen- 
dence of the ECS on the chemical potential we will fo- 
cus on the cross-section of the ECS with a (110) plane 
through the origin. This cross-section includes the com- 
plete information from all four calculated surfaces, be- 
cause they all possess surface normals within this plane. 
The ECS is shown for three different chemical environ- 
ments in Fig. pd]. Note that in a Ga-rich environment the 
(lll)(\/l9 x \7l9) reconstruction would be energetically 
more favorable than the (TTT)(2x2) Ga-adatom recon- 
struction used for the construction of the ECS at this 
chemical potential, i.e., the experimental (111) facet ap- 
pears somewhat closer to the origin. For an As-rich envi- 
ronment we find that all four considered surface orienta- 
tions exist in thermodynamic equilibrium. Furthermore, 
the (111) surface exists within the full range of accessi- 
ble chemical potentials. This is in contrast to the result 
Chctty. and Martir£3 derived from the work of Kaxiras 
et alxA: They stated that the (111) surface has a high 
energy and thus it should not exist as a thermodynamic 
equilibrium facet. However, experimental work of Weiss 
et alia using a cylindrical shaped sample indicates that 
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between the (110) and (111) orientation all surfaces facet 
into (110) and (111) orientations. The (2 x 2) super- 
structure of the (111) surface has been observed on these 
faceted surfaces. If the (111) orientation of GaAs were 
instable, the appearance of facets other than (111) on the 
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cylindrical crystal were to be expected. 

In Fig. [ll] one can see that the ECS becomes smaller 
for As-rich environments. The As terminated reconstruc- 
tions have surface energies about 20 % smaller than those 
found in Ga-rich environments, which are mostly stoi- 
chiometric like the Ga vacancy. In contrast to the sur- 
face reconstructions found for As-rich environments no 
similar Ga terminated reconstructions are observed. An- 
other remarkable feature of the ECS is that the surface 
energies do not vary very much with the orientation. For 
Ga-rich environments they vary by about ±10%, whereas 
for As-rich environments they vary only by ±5%. 

Our calculated ECS imposes restrictions on the sur- 
face energies of other surface orientations: When it has 
been proven experimentally that a facet exists in ther- 
modynamic equilibrium, one can derive a lower and an 
upper limit for its surface energy. The limits are given by 
the surface energy of the neighboring facets on our ECS 
together with appropriate geometry factors. They follow 
from the conditions that (a) the surface energy has to 
be sufficiently small, so that the surface does not facet 
into {110}, {100}, {111}, and {III} orientations, and (b) 
that the surface energy is not so small that neighboring 
facets are cut off by this plane and thus vanish from the 
ECS. In a similar way the Wulff construction yields a 
lower limit for the surface energy of any facet that does 
not exist in thermodynamic equilibrium. 

Recently the shape of large three-dimensional InAs 
islands (diameter ~ 2000 A) grown by MOVPE on a 
GaAs(100) substrate has been observed by E. Steimetz et 
aljEi These islands are presumably relaxed, the misfit of 
the lattice constants being compensated by a dislocation 
network at the InAs-GaAs interface. Thus the facets dis- 
played on these islands should be identical to the facets 
on the ECS of InAs. In fact, the observed shapes are 
compatible with an ECS like that of GaAs shown in Fig. 
|(| with {110}, {100}, {111}, and {111} facets being 
clearly discernible. Due to the similarity between InAs 
and GaAs we take this as another confirmation .of our 
results as opposed to those of Chetty and Martin.EJ 



meV/A 2 . This is due to a different splitting of the slab 
energy into contributions from the (111) and (111) sur- 
faces, as Chetty and Martin's and our sum of the (111) 
and (111) surface energies are essentially equal. Obtain- 
ing high surface energies for the (111) surfaces Chetty 
and Martin have to conclude that the (111) facet should 
be unfavorable and not exist in thermodynamic equilib- 
rium. In contrast our surface energies for the (111) sur- 
face are lower and therefore we conclude that it exists 
in thermodynamic equilibrium which appears to be in 
agreement with experimental observations. ■— ■ 

As already stated by Chetty and Martintl there are 
significant differences between the results of Kaxiras et 
al. and Northrup et al. for the (111) surface: Kaxiras et 
al. find the Ga trimer structure to be energetically fa- 
vorable in Ga-rich environments, whereas we agree with 
Northrup et al. and find it energetically unfavorable. 
This is also confirmed by experiment. 

Having calculated the absolute surface energies for dif- 
ferent orientations we are in the position to construct the 
ECS of GaAs. We have to keep in mind, however, that 
it is implicitly assumed that only the (110), (100), (111) 
and (111) surfaces exist in equilibrium. For a more re- 
fined discussion of faceting further calculations also for 
higher-index surfaces would have to be performed. From 
our ECS we conclude that in As-rich environment all four 
orientations exist in thermodynamic equilibrium. For a 
given chemical potential the variation of the surface en- 
ergy with orientation is small and less than ±10%. Our 
ECS of GaAs gives indication for the ECS of InAs or 
other III-V semiconductors which show similar surface 
reconstructions . 
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V. SUMMARY AND CONCLUSION 



The GaAs surface energies of different orientations 
have been calculated consistently with one and the same 
parameters and pseudo-potentials. The surface energies 
of the (110), (100), (111) and (111) surfaces are given in 
dependence of the chemical potentials. 

For the (111) and (111) surfaces we find a laege dif- 
ference to previous results of Chetty and MartinJla They 
derived a difference of about 62 meV/A 2 between the 
surface energies of the ideal (111) and (TTT) surfaces, 
whereas we calculate a difference of about —4 me v/A 2 . 
Consequently the absolute surface energies calculated, 
by Chetty and Martin using data of Kaxiras et alxjxB 
and Northrup et al.cB contain the above difference of 66 
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